Storing unitary operators in quantum states 
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We present a scheme to store unitary operators with self-inverse generators in quantum states and 
a general circuit to retrieve them with definite success probability. The continuous variable of the 
operator is stored in a single-qubit state and the information about the kind of the operator is 
stored in classical states with finite dimension. The probability of successful retrieval is always 1/2 
irrespective of the kind of the operator, which is proved to be maximum. In case of failure, the result 
can be corrected with additional quantum states. The retrieving circuit is almost as simple as that 
which handles only the single-qubit rotations and CNOT as the basic operations. An interactive way 
to transfer quantum dynamics, that is, to distribute naturally copy-protected programs for quantum 
computers is also presented using this scheme. 

PACS number : 03.67.Lx 



Quantum computers store information in quantum 
states and process it quantum mechanically, which make 
it possible to solve certain problems much faster than 
classical computers. The quantum states containing 
information are generally in superposed or entangled 
states that have no classical analog, and the quantum- 
mechanical processing is realized by unitary operations 
while the processing of a (reversible) classical computer 
is limited to permutation operations. Recently, schemes 
to store the processing operation itself in quantum states 
have been proposed by Preskill ^ in the context of 
the fault-tolerant quantum computation and by Vidal 
et al. from the point of view of the programmable 
gate array. In other words, quantum dynamics is stored 
in quantum states. Storing and retrieving an operation 
means that a unitary operator U, which will be applied 
to a state is encoded into a state \U), and some op- 
erator G makes the whole state evolve as G{\U) ^ \d)) — 
\U')(^{U\d)), where \U') is some residual state. The most 
distinguished difference between the program stored in 
quantum states and the one stored in classical states is 
that the former is naturally protected from copy or even 
reading. 

The previous works proposed a way to store only the 
single-qubit rotation about z axis in a quantum state and 
retrieve the operation with definite success probability. It 
might be interesting to see what other operations can be 
stored in quantum states and how they can be handled. 
Implementation of operations would be easier when one 
can store and transfer several different kinds of opera- 
tions than when one can handle only a few basic opera- 
tions and the operation of interest has to be decomposed 
first into a sequence of those basic operations. In this 
work, we generalize the scheme to store and retrieve ar- 
bitrary unitary operators satisfying some conditions and 
present a way to transfer them. 

An arbitrary unitary operator U b {d) can be written as 



trary dimension, and ^ is a real number. The number 
of distinct unitary operators is infinite because 9 is ar- 
bitrary and there are infinitely many different genera- 
tors. Therefore, it seems that infinitely many resources 
are required to store arbitrary operators The num- 
ber of different kinds of generators we need to handle is 
finite because the combinations of basic operators, such 
as the single-qubit rotations and controUed-NOT (cnot) , 
can make arbitrary unitary operators as well known [H. 
Therefore, the whole point of storing a unitary opera- 
tor is to store a continuous variable 9 in quantum states. 
Storing a real number in a digitized state requires an 
infinitely large resource, whether it is a quantum or clas- 
sical system. Since quantum system has both the digital 
and analog characteristics, however, it is possible to store 
unitary operators in finite resources as shown below, only 
if we allow the possibility of failure when retrieving. A 
finite analog system is, in principle, capable of storing 
real numbers, though there is the question of precision in 
the operation and measurement in practice. 

Consider an operator Ub{9) of an a priori known B. If 
B is not only Hermitian but also unitary, B is self-inverse 
or B^ — E where i? is a unity operator. Product oper- 
ators including Pauli operators belong to this case, 
and any Hermitian operator can be expressed as a linear 
combination of them. Then, U b (9) is decomposed as 



Ub{9) = cos{9/2)E -isin{9/2)B. 



(2) 



Here, Ub{9) is expressed as a linear combination of two 
different operators, E and S, and the information about 
9 is included in the coefficients. This expression suggests 
that Ub{9) can be stored in a single-qubit quantum state, 
say an angle state, defined by 



16*) = cos(e'/2)|0) - isin(6'/2)|l). 



(3) 



Ub{9)=cM-<0/2)B], 



(1) 



where the generator _B is a Hermitian operator of arbi- 



That is, two operators E and B are mapped onto the 
two states |0) and |1), respectively, and the coefficients 
containing 9 remain same. Note that this angle state 
contains only the information about 9, not about B at 
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all. The information about the mapping of the operators 
can be stored in additional qubits, as will be discussed 
soon. 

The operator Ub{0) is retrieved from the angle state 
by using a gate array Gb that consists of a controlled- -B 
defined by 



'E+\1){1\®B, 



(4) 



and a single-qubit Hadamard operator H on the angle 
state as shown in Fig. |^. Since B is assumed to be uni- 
tary, the controlled-i? is also a unitary operator and, 
therefore, implement able. The total dynamics of the 
joint state \6)\d) = \6) (g) \d), where the data state \d) is a 
multi-qubit state as E and B have arbitrary dimension 
in general, are described in terms oi Gb, 



GB\0)\d) = H[cos{e/2)\0)\d) -tsm{0/2)\l)B\ 

= ^[\o)UBie)\d) + \i)UB{-e)\d)]. 



(5) 



A projective measurement of the angle state in the basis 
{|0), |1)} will make the data state collapse into either a 
desired state UB{0)\d) or a wrong state UB{—0)\d) with 
the equal probability. Therefore, UB{d) stored in the 
angle state is retrieved in a probabilistic way. In case of 
failure, the correct state may be obtained by executing 
Gb once more with an additional angle state as discussed 
in Refs. If we measure the angle state after the 

execution of Gb on the joint state of a new angle state 
\29) and the wrong data state UB{—d)\d), it will give the 
desired state or a new wrong state 30)|(i) with the 

equal probability, again. This process can be repeated 
with the angle state \2™'6) [m — 2,3, •••) until we get 
the right result. The average number of the angle states 
needed for success is given by X]m=i ^^{1/2)"^ = 2. 

The simplest kind of unitary operators is the single- 
qubit rotation about a axis which is written as Eq. (|^) 
with B = Ga where CTq is a Pauli operator. In the previ- 
ous schemes, the angle of rotation was encoded in phase 
rather than amplitude like ours but they are essentially 
equivalent for single-qubit operations. One advantage of 
storing the angle information in amplitude is that the 
simple mapping in Eq. (3) and the retrieving circuit in 
Fig. 1 is applicable to any operators having self-inverse 
generators. In Ref. it is proven that the maximum 
probability of successfully retrieving the z-rotation stored 
in a single-qubit state is 1/2. It is straightforward to ex- 
tend the proof to the general unitary operators with self- 
inverse generators in our scheme (see Appendix) . There- 
fore, there exists no scheme with higher probability of 
retrieval success than ours to store general unitary oper- 
ators in a single-qubit state. 

It would be helpful to discuss what useful operators can 
be stored in the angle state and present concrete circuits 
for the gate arrays to retrieve these operators. Consider 
an n-bit coupling operator given by 



which belongs to the kind of operators of our interest 
because {(Jiz(J2z • • ■(^nzY' — E. The corresponding gate 
array Gi2 - n is illustrated in Fig. |[ This circuit con- 
sists of only n two-qubit operators and a single-qubit 
Hadamard operator because the controlled-crizCT2z • • • Cnz 
is equivalent to the n controlled-tTiz (i = 1,2, 
Therefore, as can be seen in Fig. |^, there are direct 
interactions (vertical solid lines) only between the an- 
gle qubit and each data qubit and the data qubits have 
no interactions among them, though Ji2-. „iO) itself in- 
cludes interaction among all qubits. This nice feature 
would be useful for the quantum computer using a spe- 
cial "head qubit" , which moves to mediate interactions 
between non-interacting qubits 

One of the important operators in quantum logic al- 
gebra is the controlled gate such as the Toffoli gate and 
phase-shift gate 0. For example, an n-bit phase-shift 
gate diag[l, 1, • • • , 1, e*^] is used in the quantum factoring 
algorithm. Consider the operator. 



exp 



■p: 



E 



(7) 

where a; denotes -f or — , and P^ are the projection op- 
erators (i?i±(Tiz)/2 of the i-th qubit, respectively. In the 
summation, the case of all a^'s being minus is excluded. 
This operator is equivalent to exp(— i0i7„z/2) if all the 
first {n — \) data qubits are in the state |1) and e~'^/^£^„, 
otherwise. Therefore, this is equivalent to the phase- 
shift gate up to an overall phase. Since {P^Y = -Pj^: 
P+ + P~ = E,, and P+P- = P-P+ = 0, the generator 
satisfies the self-inverse condition. Therefore, this opera- 
tor can be stored in the angle state with the correspond- 
ing controlled- -B given by an n-bit controlled-cr„2 which 
takes the angle qubit and (n — 1) data qubits as control 
bits and applies cTz to the n-th data qubit depending on 
the state of the control bits. 

We have described how a unitary operator with a self- 
inverse generator can be stored in quantum states, and 
given the corresponding gate array Gb that retrieves the 
operator from the quantum states. The angle state con- 
tains the information about 6, and the information about 
the generator B is included in the circuit Gs as the 
controlled- i?. The information about the generator is 
prerequisite for constructing the circuit. Instead of this 
impractical design, we can store the information about 
the generator somewhere else and construct a general 
circuit that interprets that information and execute the 
corresponding Gb- Number of generators to consider is 
finite for finite number of data qubits. 

One simple way of constructing the general retriev- 
ing circuit is to handle the set of the coupling operators 
Jij...k{0) defined by 



Jl2 -n{6) = exp(--z6'crizCr22 ' ' ' <7nz/2), 



(6) 



J, 



ij---k 



(0) 



• (ykzl'2), 



(8) 
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where the generator includes a subset of n spin operators 
in general. These operators include exjp{—i9f7iz/2) and 
Ji2 - n{0) in Eq. (^), and make complete set with not 
gates to produce any unitary operators such as the phase 
shift gate in Eq. (7). There are Y.Z=i nCm = (2" - 1) 
different generators of this kind. They can be stored in a 
command state consisting of log2 2" — n qubits by map- 
ping the generators onto the eigenstates of the command 
state. Since the mapping onto the eigenstates is equiva- 
lent to using classical bits, one qubit for the angle state 
and n classical bits for the command state are all that 
required to store all of these coupling operators. 

The general gate array in Fig. H is slightly modified 
from the gate array Gi2...n in Fig.H to include the com- 
mand state. For example, if Jin{0) is required to be 
stored, then the corresponding controlled- i? is the multi- 
plication of (7iz and each of which is controlled by the 
angle state. Therefore, the angle state contains 9 and the 
command state is the binary string |10 • • • 01), which indi- 
cates that only the controUed-CTiz and controlled-tT„2 are 
to be activated. To make the circuit complete, X = cr^; 
(not gate) should be added per each qubit. Since X 
is a fixed operator, one can easily include them in the 
scheme by employing one more classical bit |co) in the 
command state. For decoding, Xi is controlled by |co) as 
well as \ci) (i = 1, 2, • • • , n), which is nothing but the Tof- 
foli gate operation. Consequently, {n + 1) classical bits 
and one qubit are used to store any operators in the form 
of Jij...k{0) and Xi which can build up arbitrary unitary 
operators. The circuit is almost as simple as that which 
handles only the single-qubit rotations and CNOT. 

This scheme of storing and retrieving quantum opera- 
tions can be used to distribute naturally copy-protected 
programs for quantum computers. A programming of an 
quantum algorithm means the process of decomposing 
the unitary operation required by the algorithm into a 
sequence of basic operations. The more basic operations 
we have, the easier to program. A program can be stored 
either in classical states or quantum states. If a program 
is stored in quantum states, it can neither be copied nor 
read and only probabilistically retrievable. One way to 
distribute a quantum program is as follows. Suppose 
that Alice has her operator programed in the form of a 
sequence of the basic operators. Then, (i) Alice stores the 
first basic operator of the sequence in the angle and com- 
mand states, and sends them to Bob, (ii) Bob performs 
the gate array of Fig. ^, and tell Alice the measurement 
result - whether the operation has succeeded or not, (iii) 
Alice sends to Bob the new angle state in case of failure 
or the next operator of the sequence when succeeded, 
and (iv) they repeat (ii) and (iii) until the last operator 
of the sequence is transferred. Although there is possi- 
bility of failure in each operation transfer, it can always 
be corrected and the average number of the angle states 
necessary for successful operation is only two. This is 
an interactive distribution scheme where Alice and Bob 
have to communicate with each other during the transfer 
to guarantee the successful operation. 



In conclusion, we have shown that the unitary opera- 
tors with self-inverse generators can be stored in a quan- 
tum state and retrieved exactly by encoding the contin- 
uous variable into the probability amplitude of the state, 
at the cost of possible failure. This probabilistic feature 
is the cost we have to pay to store a real number in a 
quantum state of finite dimension, even which is impos- 
sible in classical systems. Utilizing the circuits for the 
coupling operators and storing the information about the 
generators in classical states of finite size, it is possible 
to store and retrieve arbitrary operators (including such 
ones having non-self-inverse generators), and the general 
retrieving circuit is very simple. 

This work was supported by NRL Program, electron 
spin science center, and BK21 Project. 



APPENDIX: 

Any scheme to store the quantum operation U{9) = 
exp[-|-i(6'/2)i?] of a fixed, self-inverse generator B in the 
quantum states can be described by a unitary transfor- 
mation, 



G{\Ug) ® \d)) = 



®u{e)\d) + Ji-pi\xf) 



(Al) 



where \Ue) is the program state containing the informa- 
tion about the operation and \d) is the data state. G 
transforms the total state into the sum of the product 
of some residual state |Tg) and the desired state U{9)\d) 
with success probability pg and the failed state \xg) with 
probability 1—pg. To distinguish success from failure by 
measurement, it should be satisfied that (t^ Ixg/) = for 
all d, d', 6*, and 9'. 

Suppose that \d) is an n-qubit state and iV = 2". The 
data state \d) is expanded by the computational basis 
as Id) — '^k^QCk\k), where Ck is the complex coeffi- 
cient that satisfies the normalization condition. Then, 
Eq. ( |Al| ) is rewritten as 



G{^Ck\Ug)\k) 



)U{9)\k) 



(A2) 



where product sign (g) between the program and data 
states is omitted for simplicity. This implies that pi and 



\Tg) do not depend on \d) because RHS's of Eqs. ( [Al| ) 



and (|AJ) must be same for all \d). Now, pg and \Tg) will 
be denoted by pg and |re), respectively. 

The one-qubit program state can be also expanded 
by \Ue) = a{9)\0) + P{9)\tt_), where 0(6*) and P{9) are 
complex function satisfying the normalization condition 
{Ug\Ug) — 1, and |0) and |7r) are the program states cor- 
responding to the operations E and i?, respectively. For 
any working scheme G, there always exist the program 
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states for = and 9 = tt. These states correspond to 
E and B, respectively, because U{9) = exp[+i{9/2)B] = 
cos{9/2)E + ism{9/2)B, Ue = E ior 9 = 0, and Ug = B 
for 9 — TT (up to an overall phase). Therefore, we can 
always expand \Ue) by the linear combination of |0) and 
Itt), which are not necessarily orthonormal to each other. 
Again, Eq. (Al) is expressed as 



Giamm) + m\K)\d)] = 

a{9) [^\To)E\d) + v/l-po|Xo> 

+ m \^\n)B\d) + v/T^lx?) 



(A3) 



for all Id). Therefore, all \Tg), \to), and |ri) are same with 
|r) not depending on and 



mU{9) = a{9)^E + mVP^B. 



(A4) 



This means that a{9) = ^/pgJpQ cos{9 / 2) and /3{9) = 
\/pe/pi sm{9/2). The remaining part of the proof is same 
with that of Ref. |^ . The maximum probability of success 
is 1/2, which is achieved by our scheme. 
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FIG. 2. The gate array G\2 -n for the n-bit couphng oper- 
ator Ji2.- n(S). Zi represents for the i-th data qubit. 
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FIG. 3. The general gate array handling a set of basic op- 
erators, Jij.-.kifi) and Xi. 
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FIG. 1. The general gate array Gb that retrieves Ub(Q) of 
a fixed generator B from the angle state \9). H is a Hadamard 
operator and M represents a projective measurement in the 
basis {|0), |1)}. The short diagonal line on the bottom line 
represents that the data state |d) consists of several qubits. 
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